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Abstract. In this review I discuss the characteristics and the formation of all classes 
[t"! I of millisecond pulsars (MSPs). The main focus is on the stellar astrophysics of X- 

HH , ray binaries leading to the production of fully recycled MSPs with white dwarf (WD) 

^"i • or substellar semi-degenerate companions. Depending on the nature of the compan- 

' ion star MSPs are believed to form from either low-mass X-ray binaries (LMXBs) or 

Qh. intermediate-mass X-ray binaries (IMXBs). For each of these two classes of X-ray bi- 

naries the evolutionary status of the donor star - or equivalently, the orbital period - at 
^ I the onset of the Roche-lobe overflow (RLO) is the determining factor for the outcome 

. of the mass-transfer phase and thus the nature of the MSP formed. Furthermore, the 

' formation of binary MSPs is discussed in context of the PP-diagram, as well as new 

interpretations of the Corbet diagram. Finally, I present new models of Case A RLO of 
IMXBs in order to reproduce the two solar mass pulsar PSR J1614-2230. 
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1. Introduction 



MSPs are known to be important sources of research in many areas of fundamental 
physics. Equally important, however, binary MSPs represent the end point of stellar 
evolution, and their observed orbital and stellar properties are fossil records of their 
kJi • evolutionary history. Thus one can also use binary pulsar systems as key probes of 

', Stellar astrophysics. It is well established that the neutron star in binary MSP sys- 

tems forms first, descending from the initially more massive of the two binary stellar 
components. The neutron star is subsequently spun-up to a high spin freq uency via 
accretion of mass and angular momentum once the secondary star evolves (Alpar et alj 
IT982 : ,Radhakrishnan & Srinivasan,1982. : .Bhattacharya & van den Heuvel,1991,). Dur- 



ing this recycling ph ase the system is observable as an X-ray binary (e.g. Nagase|[l989l : 
iBildsten et al. 199 7|) and towards the end of this phase as an X-ray millisecond pul- 
sar ( Wijnands & van derKlisllT998l : I Archibald et aDl2009h . Although this formation 



scenario is now commonly accepted, many aspects of the mass-transfer process and 
the accretion physics (e.g. details of non-conservative evolution and the accretion ef- 
ficiency) are still not well understood. Some of these issues will be addressed here. 
Given the limited number of pages in this conference proceedings review the discus- 
sion will mainly be of a qualit a tive c h aracter. For a more deta i led re view, see e.g. 



Bhattacharya & van den Heuveil (ll99lh : lTauris & van den Heuvell (l2006l) . 
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2. Mass transfer in X-ray binaries and the nature of the donor star 



Consider a close interacting binary system which consists of a non-degenerate donor 
stai" and a compact object, in our case a neutron star. If the orbital separation is small 
enough the (evolved) non-degenerate star fills its inner common equipotential surface 
(Roche-lobe) and becomes a donor star for a subsequent epoch of mass transfer toward 
the, now, accreting neutron star. In this phase the system is observed as an X-ray bi- 
nary. When the donor star fills its Roche-lobe it is perturbed by removal of mass and 
it falls out of hydrostatic and thermal equilibrium. In the process of re-establishing 
equilibrium the star will either grow or shrink - depending on the properties of its en- 
velope layers - first on a dynamical (adiabatic) timescale and subsequently on a slower 
thermal timescale. However, any exchange and loss of mass in such an X-ray binary 
system will also lead to alterations of the orbital dynamics via modifications in the 
orbital angular momentum, and hence changes in the size of the critical Roche-lobe 
radius of the donor star. The stability of the mass-transfer process therefore depends 
on how these two radii evolve (i.e. the radius of the star and the Roche-lobe radius). 
The various possible modes of mass exchange and loss include, for example, direct 
fast wind mass loss, Roche-lobe overflow (with or without isotropic re-emiss ion) and 
common envelope evolution (e.g. lvan den Heuvellll994l : ISoberman et al. I ll 9971. and ref- 
erences therein). The RLO mass transfer can be initiated while the donor star is still 
on the main sequence (Case A RLO), during hydrogen shell burning (Case B RLO) or 
during helium shell burning (Case C RLO) - see Fig. [T] The corresponding evolution- 
ary timescales for these different cases will in general proceed on a nuclear, thermal or 
dynamical timescale, respectively, or a combination thereof. This timescale is impor- 
tant for the extent to which an old neutron star can be recycled (i.e. with respect to its 
final spin period and B -field). 
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Figure 1 . Stellar radius as a function of age f or a 5 Mp star. Also shown is th e 
definition of RLO Cases A, B and C. Figure from lTauris & van den Heuvell ( l2006l) . 
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The dynamical evolution of a binary system can be found by solving for the 
changes in the orbital separation, a. The orbital angular momentum of a circular bi- 
nary system is given by: /orb = i^0.a^, where /u is the reduced mass and the orbital 

angular velocity is: Q = yjGM/a^. A simple logarithmic differentiation of the orbital 
angular momentum equation yields the rate of change in orbital separation: 

where the two stellar masses are given by M\ and M2, and the total change in orbital an- 
gular momentum per unit time is given by: 7orb - jgm+jmb+j\?,+Jm\- Thcsc four tcrms 
represent gravitational wave radiation, magnetic brakin g, spin-orbit couplings and mass 
loss, respectively (e.g. 'Taur is & van den Heuvelll2006l and references therein). 

When modelling the evolution of an X-ray binary one should take into account 
a number of issues related to: the stellar evolution, the stability of the mass-transfer 
process, the ejection of matter from the system and the accretion of material onto the 
neutron star - see illustration in Fig. |2] Ideally, all these calculations should be per- 
formed self-consistently. The major uncertainties here are related to the amount and 
the mode of the specific orbital angular momentum of ejected matter - and for close 
systems, also the treatment of the spin-orbit couplings. 
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Figure 2. 



Illustration of some of the many issues to consider in X-ray binary models. 
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Figure 3. The distribution of the ~ 100 binary radio pulsars with measured values 
in the PP-diagram. The symbols to the right indicate the nature of the companion 
star. Pulsars with a He WD companion marked by a ring are discusse d in Section|6] 
Data taken from the ATNF Pulsar Catalogue in May 201 1 dManchester et alj|2005i 
|http://www.atnf.csiro.au/research/pulsar/psrcatj . 



3. The population of binary millisecond pulsars 

Fig. [3] shows the location of the known binary pulsars in the PP-diagram. The fully 
recycled MSPs (P < 10 ms) are dominated by having mainly He WD companions, 
although also ultra light (substellar) semi-degenerate companions are seen - often in 
eclipsing "black widow-like systems" - as well as a few systems with the more massive 
CO WD companions. The mildly recycled MSPs (10 ms < P < 100 ms) are dominated 
by CO WD (or ONeMg WD) and neutron star companions. As we shall see, this rela- 
tively slow spin rate is expected from an evolutionary point of view as a consequence 
of the rapid mass-transfer phase from a relatively massive donor star. The pulsars with 
similar spin periods and He WD companions apparently also had a limited recycling 
which may hint their origin - see Section |6] The double neutron star systems descend 
from high-mass X-ray binaries (HMXBs) and are not disc ussed in this review (see e.g. 
Voss & Taurisll2003l : iDewi & Polsll2003l : iDewi et al.ll200^ . and references therein). 



MSPs are not only characterized by a rapid spin. All of them also possess a low surface 
magnetic flux density, B which is typically of the order 10^ G, or some 3-5 orders of 
magnitudes less than the B -fields of ordinary, non-recycled pulsars (see discussion in 
Section |7]). Hence, the recycled pulsars do not suffer as much from loss of rotational 
energy due to emission of magnetodipole waves. For this reason the MSPs have small 
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period derivatives P < 10"'^ and hence they are able to maintain the production of 
radio waves keeping them observable for billions of yeai^s. 

Binary pulsars in globular clusters are, in general, not suitable as tracers of their stellar 
evolution history because of the frequent encount ers and exchanges of companion stars 
in the dense environments (Ransom et al.l l2005l) . Hence, this review only considers 
binary pulsars formed in the Galactic disk. 

When discussing the origin of MSPs with various companion stars it makes sense 
to plot the binary orbital period as a function of the companion star mass - see Fig. |4] 
This figure is essential for understanding the progenitor systems. 
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Figure 4. The distribution of 162 binary radio pulsars in the companion mass- 
orbital period plane. Data taken from the ATNF Pulsar Catalogue in May 201 1 . The 
five roads to form an MSP are shown: Two roads from LMXBs (Cases A and B) and 
three roads from IMXBs (Cases A, B and C), see also Fig. Hand Tabled HMXBs 
leave behind the double neutron star systems and only partially recycle the first bom 
neutron star (see their location in Fig. |3). The boundaries are not very strict and keep 
in mind that most of the companion masses have large error bars. 



4. MSPs with He WD or substellar dwarf companions - evolution from LMXBs 

To the left in Fig. H] one sees the MSPs with either ultra light companions or He WD 
companions. Having low-mass companions these systems thus descend from LMXBs. 
MSPs with Porb < 1 day originate from LMXBs in very tight orbits where the donor 
star akeady initiated RLO while it was still on the main sequence (thus Case A RLO), 
whereas the MSPs with Pq* > 1 day originate from wider orbit LMXBs where the 
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donor star did not fill its Roche-lobe until it had ev olved and expanded to become a 
(sub)giant, i.e. Case B RLO. It has been shown by 'Pylyser & Savonij^ (I1988L Il989l) 
that a critical orbital bifurcation period {Pbn) exists at the onset of the RLO, separating 
the formation of converging LMXBs from diverging LMXBs. The theoretical estimated 
value of Pbif is ~ 1 day, but depends strongly on the treatment of tidal interactions and 
the assumed strength of magnetic braking which dra ins the system of orbital angular 
momentum dvan der Sluvs et al.ll2005l : iMa & Lill2009h . 



4.1. Close LMXB systems - formation of black-widow pulsars 

In observed MSP systems with forb ^ 1 day the mass transfer was driven by loss of or- 
bital angular momentum due to magnetic braking and (for more narrow systems) emis- 
sion of gravitational waves. The donor star becomes degenerate, or semi-degenerate, 
as the ultra-compact binary evolves. An additional mass-transfer phase is initiated 
in those final NS-i-WD systems where emission of gravitational waves drives the two 
stellar components together. The outcome is ultra-compact systems (e.g. X-ray tran- 
sient accre ting MSPs) with orbital periods less than one hour and M2 - 10"^ M© 



dPeloye & Bildsten 2003). Furthe rmore, irradiation effects, first by the X-ray flux 
fr om the accreting neutron star (jPodsiadlowskil Il99lb and later by the pulsar wind 

fivani & Brookshaw"l992f), as well as tidal dissipation of energy in the envelope 
pplegate & Shaham 1994), may cause the companion star to be thermally bloated 
and evaporate. In some cases this combi ned evolution eventually leads to the forma- 
tion of pulsar pla nets (PSR B 1257-H12, Wolszczan & Fraillll992h and sohtary MSPs 
(PSR B 1937-1-21. [Backer et al.| [l982). Many radio pulsars in narrow systems display 
eclipses during their orbital motion and in the case of PSR J20 5 1-0827 one can eve n 
measure the effects of gravitational quadrupole moment changes ( Lazaridis et al.l20Tlh . 
In those binaries where the companion leaves behind a low-mass He WD (< 0.20 Mq) 
the cooling age determination of the WD (and thus the age of the MSP) is complicated 
by their thick residual hydrogen envelopes. Pycno-nuclear burning at the bottom of 
these envelopes can keep the low-mass He-WDs warm for ~ 10^ yr (see discussion in 
Ivan Kerkwijk et al.ll200a) . 



4.2. Wide LMXB systems - formation of classic MSPs with He WD companions 

In wide-orbit LMXBs the donor star did not fill its Roche-lobe until it moved up the 
red giant branch (RGB). For low-mass stars (< 2.3 M©) on the RGB there is a well- 
known relationship between the mass of the degenerate helium core and the radius 
of the giant st ar - almost entirely independent of the mass present in the hydrogen- 
rich envelope (iRefsdal & Weigert 197lUWebbink et aLlll983 ). This relationship is very 
important for the formation of wide-orbit MSPs since it results in a unique relation- 
ship between orbital period and white dwarf m ass (.Savoniiel [l987l : iJoss et al.lll987l : 
Rappaport et all [19951 : iTauris & Savonijel Il999l) . The companions here are He WDs 
with masses 0.2 < Mwd /Mq < 0.46 - unless the system was initially so wide that the 
donor ascended the asymptotic giant branch (AGB) before initiating RLO and thus left 
behind a CO WD companion (as observed only in PSR B0820-I-02 which has an or- 
bital period of 1232 days). The correlation between Porb and Mwd is difficult to verify 
observationally since very few MSPs have accurately measured masses of their com- 
panion. From the observed mass function one can only estimate the WD mass which 
depends on the orbital inclination angle of the system as well as the neutron star mass 
- both of which are unknown. The result of assuming a fixed neutron star mass which 
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is the same for all systems and plotting the orbital periods as a function of the median 
He WD masses (i.e. assuming an inclination angle, / = 60° for all the systems) can be 
seen in Fig. [5] At first sight this correlation is not obvious under these assumptions - not 
even when considering that the correlation depends on the chemical composition of the 
donor star as well as the treatment of mixing. However, van Kerkwij k et all (l2005h have 
shown that the relation is indeed quite good if one only considers those 7-8 systems for 
which the WD mass has been estimated fairly accurately. On the other hand, there 
might be a systematic deviation from the correlation for pulsars with PqA > 100 days 
(ffauris 1996 ; Tauris & Savonijelil999; Stairs et al. .2005) . although this claim is based 
on very small number statistics. Hopefully, future mass determinations of some of these 
wide binaries can help to settle this issue. 

Finally, it should be mentioned that these wide-orbit binary MSPs also have an- 
other fossil of the mass-transfer phase: the correlation between orbital eccentricity and 
orbital period. These so-called residual eccentricitie s are, in general, very small (typi- 
cally between 10"^ - 10"^). The correlation (P hinneyll 19921) is related to tides and arises 
because density fluctuations in the (convective) envelope increase with more evolved 
donor stars, which have wider orbits, thus preventing perfect circularization. 




0.5 



Mass of white dwarf, (Mq) 



Figure 5. Orbital period as a function of WD mass for MSPs with He WD com- 
panions. The error bars of the WD masses represent the unknown orbital inclination 
angle - the left end corresponds to an inclination angle of 90° and the right end 
marks the 90% probability limit. The additional assumption is that all neutron stars 
are assumed to have a mass of 1.60Mq (for simplicity). This mass is larger than 
the typical ~ 1.35 Mq found for pulsars in double neutron star systems. However, 
fully recycled MSPs are expected to have accreted signifi cantly more material. Th e 
theoretical (Mwd, /'orb)-correlation is shown as a blue line (iTauris & Savoniiell999h . 
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5. MSPs with CO WD companions - evolution from IMXBs 

CO WDs (and ONeMg WDs) are substantially more massive than He WDs. Therefore, 
MSPs with these more massive WDs must originate from binaries with more massive 
WD progenitor stars - i.e. IMXBs which have don or mas ses of typically 3-6 Mq. 
Observationally, these systems were first identified by lCamilo. (1996) as an independent 
class of MSPs. There seems to be three roads to produce such systems (see below). 

5.1. Wide IMXB systems - Case C RLO and common envelope evolution 

It was demonstrated nicely by Ivan den Heuvell ( 1994 ) that binary systems with a mildly 



recycled MSP and a CO WD companion could form from an IMXB with a donor star 
on the AGB, leading to common envelope (CE) and spiral-in evolution once the RLO 
sets in. This model is in particular the favoured scenario for the formation of very tight 
binary MSPs with orbital periods Porb < 3 days. The reason for this is the ability of the 
spiral-in phase to reduce the orbital angular momentum, and thus the orbital period, by 
a huge amount (e.g. Dewi & Tauris 2000. and references therein). One point to notice 
from this model is that the CE phase is expected to be extremely short. The in-spiral 
proceeds on a dynamical timescale (a few orbital periods) followed by an envelope ejec- 
tion phase that may take up to 10^^ years. However, at least ~ 0.01 Mq is needed to spin 
up a pulsar to about 10 ms (see Section IT]) and this requires some 0.4 Myr of efficient 
accretion at the Eddington limit (a few times 10~*^ Moyr"^). It is therefore believed 
that the recycling of the MSP may actually take place from either wind accretion from 
the naked post-CE core, or from Case BB RLO if the naked helium star expands and 
fills its Roche-lobe (again) leading to additional recycling of the neutron star. 

5.2. Hertzsprung gap IMXB systems - early Case B RLO and isotropic re-emission 

An alternative way of producing MSPs with CO WD companions was demonstrated by 
iTauris etaTI (12000'). They considered donor stars which had just left the main sequence 
(i.e. Hertzsprung gap, or sub-giant stars) and applied the so-called isotropic re-emission 
model of an IMXB. In this model ( Bhattacharva & van den Heuvel 119911) . the far ma- 



jority of the matter - being transfered at a highly super-Eddington rate - is ejected (e.g. 
in a jet) with the specific orbital angular momentum of the accreting neutron star. This 
model can stabilize the RLO and prevent it from becoming dynamically unstable. The 
typical donor star masses in this scenario are 3 - 6 Mq and they leave behind a CO WD. 
Hence, the outcome of this formation scenario is mildly recycled MSPs with CO WD 
companions and orbital periods between 3 and 50 days. 

5.3. Close IMXB systems - Case A RLO and a long-lasting, stable mass transfer 

The formation of MSPs from close-orbit IMXBs with main sequence donors has been 



la in seq uenc 

studied in detail by iPodsiadlowski et all (l2002l ): lLin et all (l20Tl]) : lT^uris et all (1201 lal) : 



Tauris & Langej ( 2011 ). The outcome of these systems is mainly MSPs with CO WD 
companions (but He WD companions are also formed if the donor mass is at the low 
end or the initial binary is very close). The MSPs formed through this channel are 
expected to be fully recycled (i.e. with spin periods of a few ms) since the last part of 
the mass-transfer phase lasts about 10' yr in this case (see Fig. O. This is opposite to 
the preferentially mildly recycled MSPs formed by IMXBs via early Case B or Case C 
RLO (described above) where the duration of the total accretion phase is only of the 
order of 10^ yr, or even less for Case C RLO which always evolves through a CE. 
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Figure 6. The evolution of an IMXB system undergoing Case A RLO. These 
models were calculated in order to reproduce the record high -mass (1.97 M©) MSP 
J1614-2230. The mass transfer takes place when the donor star in the HR-diagram 
(upper left panel) is on the track between the open squares, including the loop down 
to log Teff - 3.7. Three phases of mass transfer (Al, A2 and AB) are identified 
within this track - see two bottom panels which show the mass-transfer rate as a 
function of age and donor star mass, respectively. The neutron star is spun up to an 
MSP during phase AB where it accretes the far majority of its accumulated material 
(upper right panel). We found that this n eutron star was bom with an initial mass of 
1 .7 ± 0.15 Mq - see lTauris et al.l (1201 la') for a detailed description. 



5.3.1. PSR J1614-2230 

To shown an example of an IMXB system of Case A RLO we have shown in Fig. |6] 
one possible formation path for this recently discovered 1 .97 M© MSP ("Pemores t et all 
12010) . It is interesting to notice that this MSP is the only system with a CO WD 
companion which is fully recycled {P - 3.15 ms). However, this full recycling of the 
pulsar is indeed expected if it evolved from an IMXB of Case A RLO, see above. 



5.4. AIC - accretion induced collapse (a half road to form a millisecond pulsar?) 

It has been suggested that a neutron star can form from the implosion of a massive 
ONeMg WD if such a WD is accreting mater i al (within a certain rate) from a compan- 
ion star (e.g. lMichellll987l : ICanal et al. 119901 : iNomoto & Konddll99lL and references 
therein). This companion can either be a main sequence star (super-soft X-ray source), 
a low-mass giant (novae-like system) or a helium star. The AIC route to form MSPs 
has two main advantages: 1) it may help to explain the postulated birthrate problem 



10 



Thomas M. Tauris 



dKulkami & Narayanlll988l) between the number of LMXB progenitor systems and the 
observed number of MSPs, and 2) it can retain pulsars in globular clusters due to the 
small momentum kick expected to be associat ed with the irnplosion . (For details of sim- 
ulations of the implosion mechanism, see e.g. IDessart et al. liooi). However, it seems 
difficult to predict the pulsar spin rate, as well as the surface magnetic field strength, 
associated with a pulsar formed via AIC. On the other hand, one can argue that the 
newborn neutron star formed via AIC might at a later stage - once its donor star has 
recovered from the dynamical effects of sudden mass loss (caused by the released grav- 
itational binding energy in the transition from an accreting WD to a more compact 
neutron star) - begin to accrete further material from its companion star, which should 
then resemble the conditions under which MS Ps are formed via the conventional chan- 
nels. A recent population synthesis study by iHurley et al.l (12010) concludes that one 
cannot ignore the AIC route. However, there are still uncertainties involved in such 
studies and in particular the applied conditions for making the ONeMg WD mass grow 
sufficiently. A weakness in advocating the AIC formation channel is thus exactly the 
point that one cannot easily distinguish the result of this path from the standard scenario 
(see, however. Section |6]for a novel speculative hypothesis). It would be very interest- 
ing though, if observations would yield either a very slowly spinning (few hundred ms) 
pulsar associated with a very low B-field, or a high B-field MSP with a very rapid spin. 
(However, the latter kind of pulsar would be very unlikely to be detected given that its 
very strong magnetic dipole radiation would slow down its spin rate within a few Myr). 



Table 1. Progenitor systems leading to the formation of MSPs - see also Fig. |7] 
All the characteristic values for both the X-binaries and the MSP systems are only 
rough indications and depend on effects which are poorly known - for example, 
concerning the strength of magnetic braking and other spin-orbit couplings, as well 
as the CE and spiral-in evolution. AICs are not included in this table. 
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Figure 7. Formation of MSPs. From left to right: LMXBs (Cases A and B), IMXBs (Cases A, early B and C (CE)) and possibly AIC? 
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6. The Corbet diagram revisited - a new puzzle? (related to AIC?) 

Plotting all binary pulsars in the Corbet diagram does not reveal much interesting in- 
formation - it is a big mess ! However, if one only considers the population of pulsars 
with He WD companions an interesting pattern is noticed, as shown in Fig.[S] Region I 
shows that MSPs can be fully recycled over a spread of 3 orders of magnitude in final 
orbital period. (Note, the well-known orbital period gap in region I, currently between 
23 and 53 days). If the orbital period exceeds ~ 200 days the pulsars are only partially 
recycled as noticed from their slow spin periods (region III). This is expected since the 
mass-transfer phase is relatively short in those w ide LMXBs where the d onor star is 
highly evolved by the time it fills its Roche-lobe ( Tauris & Savonijdll999l) . This may 



explain why no pulsars are seen in region IV. In region 11 one sees a sub-population 
of systems with orbital periods between 1 and 20 days, all of which are only mildly 
recycled MSPs with spin periods between 20 and 100 ms. Where do these systems 
come from? Obviously from progenitor systems where the mass transfer was limited. 
Whereas the populations in regions I and III are produced from LMXBs with very dif- 
ferent initial periods the situation might be different for region II. Could these systems 
then perhaps originate from AICs where the subsequent spin up of the newborn neutron 
star only resulted in a mild spin-up? (because of the limited material remaining in the 
donor star envelope following the mass transfer to the ONeMg WD before the implo- 
sion). The narrow range of orbital periods for these pulsars could then reflect the fine- 
tuned interval of allowed mass-transfer rates needed for the progenitor ONeMg WD to 
accrete and grow in mass to the Chandrasekhar limit before its implosion. 
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Figure 8. The distribution of the 47 binary radio pulsars with a He WD compan- 
ion in the Corbet diagram. The plot reveals four regions marked by I, II, III and IV. 
These regions may be understood from an evolutionary point of view - see text.. The 
pulsars in region II are also marked by a ring in the f P-diagram in Fig. [3] 
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7. Spinning up old neutron stars to millisecond pulsars 

The physics of the accretion disk-magnetosphere interaction is still not known in de- 
tail. The interplay between the neutron star magnetic field and the conducting plasma 
in the disk is a rather c omplex process. For details of the accretion physics I refer to 
Pringle &Rees (1972) Lamb et al l(fl973h:IShapiro & Teukolsk v ( 1983): G hosh & LambI 



(120061) and references therein. When the conditions for accretion are optimal the spin- 
up torque, N = dJ/dt, acting at the edge of the co-rotating magnetosphere (roughly 
at the distance of the Alfven radius, r^) is easily estimated from the specific orbital 



angular momentum at this location {j = \ fxv\ 



d 



■\JGMrA) and is roughly: 



N ^ -{m)~ My/GMrAf 



(2) 



where M is the mass of the neutron star, M is its accretion rate and / is a factor of 
order unity depending on the accretion flow. Observations of the spin evolution of 
accreting neutro n stars at a given in stant show rather stochastic variations, even with 
torque reversals (iBildsten et al.lll997h . However, the long-term spin rate will eventually 
tend toward the equilibrium spin period, Peq- The location of the associated so-called 
spin-up line for the rejuvenated pulsar in the PP-diagram can be found by considering 
the equilibrium configuration when the angular velocity of the neutron star is equal to 
the Keplerian angular velocity of the magnetosphere, at roughly the Alfven surface, 
where the accreted matter enters the magnetosphere, i.e. = Qeq - ^^k(''a) or: 



eq 



(3) 



where coc is the so-called critical fastness parameter which is a measure of where the 
accret ion torque vanishes (depending on the dynamical importance of the pulsar spin 
rate). iTauris et all ( 201 Ibh recently demonstrated that one can obtain a convenient ex- 
pression for the amount of accreted mass needed to spin up a pulsar: 



AMeq - 0.22 Mo 



(M/Mq) 

n4/3 



1/3 



(4) 



where M is the mass of the recycled pulsar (after spin-up) and the equilibrium spin 
period is given in ms. In Fig. |9]we show how much mass is needed to spin up a pulsar 
to a given spin period. For example, it requires 0.003, 0.01 and 0.10 Mq to spin up a 
pulsar to a spin period of 25, 10 and 2 ms, respectively. The value of AMgq should be 
regarded as a lower limit to the actual amount of material transfered since the accretion 
efficiency is less than unity - even at sub-Eddington levels, see Section [8^ 

It is a fact that the B -field of the neutron star decreases during the recycling pro- 
cess. However, it is still not understood if this is caused by spin-down induced flux 
expulsion of the core proton fluxoids (ISrinivasan et al.lll990l) . or if the B -field is con- 
fined to the crustal regions and decays due to diffusion and Ohmic dissipation, as a 
result of a decreased electrical conductivity when heating effects set in from nuclea r 
burning of the accreted material (iGeppert & Urpinlll994l: iKonar & Bhattacharyalll997 ), 
or if the field decay is simply caused by a diamagnetic screening by the accreted plasma 



14 



Thomas M. Tauris 




1 10 100 

Equilibrium spin period, P^^ (ms) 

Figure 9. The amount of mass needed to spin up a pulsar as a function of its 
equilibrium spin period - see text for discussion. The curves correspond to different 
neutron star masses between 1 .0 - 2.0 Mq (in steps of 0.2 Mq) increasing upwards. 



- see iBhattacharyal ( 2002 ) and references therein. The decay of the crustal B -field is 
found by solving the induction equation: 



— = V X (- X V X B) + V X (V X B) 

ot An cr 



(5) 



where the electrical conductivity, cr depends on the temperature of the crust (as well as 
the local density and lattice impurities). 



8. Masses of millisecond pulsars 

In order to weigh a pulsar it must be a member of a binary system. The most pre- 
cisely measured masses of pulsars are obtained via general relativistic effects. The 
related post-Keplerian parameters include periastron advance, redshift/time dilation, 
orbital period derivative and Shapiro delay (e.g. Will 2009). Shapiro delays of radio 
signals from pulsars (IStairs et al.lll99 8) have the advantage of being measurable also in 
low eccentricity systems if the orbital inclination is high (i.e. close to 90°). This method 
yields the opportunity to weigh both compact stars accurately. So far, s uch measure- 
ment s have revealed masses of recycled n eutron stars between 1.17 Mq dJanssen et al. 



l2008h and 1.97 Mq ( Demorest et al.ll2010l) . This wide spread of recycled pulsar masses 



reveals both the various degree of recycling and a spread in the birth masses of these 
neutron stars - see extended discussion in [Tauris et al.l(l201 lal) . 

8.1. An anti-correlation between pulsar mass and orbital period? 

For pulsars with He WD companions it has been argued that an anti-correlation between 
pulsar mass and orbital period would be expected as a simple consequence of the inter- 
play between mass-transfer rate (and thus accretion rate), orbital period and the evolu- 
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tionar y status of the LMXB donor star at the onset of the RLO (see iTauris & Savonije 



19991 . for details). So far, this anti-correlation has not been confirmed from the few 



systems with measured pulsar masses. Keep also in mind that this anti-correlation is 
probably blurred by the fact that neutron stars are bom with different masses. 

8.2. The accretion efficiency 

Tauris&Savoniiel(ll999h noticed from a comparison between observational constraints 



on pulsar masses and their calculated LMXB models that the accretion efficiency is 
sometimes less than 50% - even for those LMXBs where the neutron star accretes at 
sub-Eddington rates. The conclusion is that a significant amount of matter may be 
ejected from the pulsar magnetosphere due to magnetodipole wave press ure, the Gunn- 
Ostriker mechanism or the propeller effect ('lUarionov & Sunyaev Il975h . Furthermore, 



accretion disk instabilities (Pringle 1981 ; van Paradijs 1996 ) are also responsible for 
ejecting part of the transfered material. Hopefully, future measurements of neutron star 
masses can help us to constrain the accretion efficiency better. 



9. Summary 

When trying to determine the progenitor system of a given observed binary MSP there 
are a number of important parameters to consider: the pulsar spin period, its period 
derivative, the orbital period, the eccentricity and the masses of both the companion 
star and the pulsar itself. Given these properties, it should be possible to place the MSP 
among the five (and a half) formation channels summarized in this paper. (The origin of 
the double neutron star systems is not considered in this review, but has an equally inter- 
esting story). Of course, there are constantly discovered new oddballs trying to destroy 
our general picture of MSP formation, calling for a continuous revision of possibil- 
ities. For example, the discovery of the binary MSP J1903-I-0327 (Freire et al. 20111) 
which requires a triple system origin to explain its current peculiar characteristics (fully 
recycled, high eccentricity, wide orbit, 1 M© main sequence companion star). 

The most important parameters for forecasting the outcome of an X-ray binary 
system are the mass ratio of the stellar components and the evolutionary status of the 
donor star at the onset of RLO. The nature of the observed MSP companion star tells 
us if the progenitor systems was an LMXB or an IMXB. 

The MSPs which are only mildly recycled must have had a limited mass-transfer 
phase: either because their companion star at this stage was relative massive and/or 
evolved enough to have a convective envelope, or it might be connected to an AIC 
event. The Corbet diagram for these binary pulsars with He WD companions seems to 
include important informations about the origin of these mildly recycled MSPs. 

It is of great theoretical desire to have measured as many pulsar and companion 
masses as possible. These values are important for understanding neutron star forma- 
tion, binary stellar evolution and accretion processes. 

It was my intention to describe the beautiful link between observations of MSPs 
and binary stellar astrophysics. If the reader finds this paper somewhat boring it is 
certainly not the fault of the topic. 
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